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ABSTRACT 


Investigators in several countries have been studying a two-frequency 
microwave radar technique for the remote-sensing of directional ocean wave 
spectra and surface currents. This technique is conceptually attractive because 
its operational physical principle involves a spacial electromagnetic scattering 
resonar e with a single, but selectable, long gravity wave. Multiplexing of 
signals having different spacing of the two transmitted frequencies allows 
measurements of the entire long wave ocean spectrun to be carried out. The 
two- frequency scatterometer signal/background ratio is, using conventional 
processing, only adequate for measurements made at or relatively near the 
dominant wave wavenumber when sufficient temporal averaging of spectra is 
employed . 

A new scatterometer has been developed and experimentally tested which is 
capable of making measurements having much larger signal/background values than 
were previously possible. This new instrunent couples in a hybrid fashion the 
resonance technique with coherent, frequency-agility radar capabilities. This 
scatterometer is presently configured for supporting a program of 
surface-current measurements. Straightforward system modifications will allow 
directional measurements of ocean wave slope spectra to be carried out as well 
as current measurements. 


1. INTRODUCTION 

The operation of dual- frequency scatterometers in detecting microwave 
backscatter from the ocean surface is now well understood. Such scatterometers 
transmit two microwave signals separated in frequency by some Af which is in the 
Megahertz range. Return signals due to each of the transmitted signals are 
separately received and then multiplied together. The power spectrun of the 
resultant output appears as shown schematically in Fig. 1. It consists of the 
sum of a broad background spectrun and a sharp "Ak-line". The background 
spectrum is the result of a convolution of the Doppler spectra of the two 
received microwave signal* and consequently has a width of a few tens to a few 
hundred Hertz. The Ak-line, on the other hand, is the result of a resonance 
between the beat pattern of the two electromagnetic signals and the modulation 
pattern of the short ocean surface waves (centimetric "Bragg waves") responsible 
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«■ for the backscatter. It has a width of a few tenths to a few Hertz depending on 

the relative speed of the microwave antenna and the ocean surface. (Plant, 

1977; Alpers and Hasselman, 1978; Jackson, 1931 ). 

* The frequency of the center of the £k-line is exactly equal to the 

< frequency of the ocean surface wave which travels along the horizontal radar 
4 look direction and whose wavenumber, K w , is 2tk cos 0 where Ak is the wavenumber 

separation of the transmitted microwave signals and 9 is the grazing angle 
► (Plant, 1977). If a current exists on the surface, the wave of wavenumber K 

will be advected by the current and its frequency in the antenna's frame of 
reference will change. By this means su-face currents along the horizontal 
radar look direction can be measured using dual-frequency scatterometers 
(Schuler, 1978; Alpers et. al., 1980). By varying Ak, vertical profiles of 
horizontal currents along the radar look direction can also be determined 
(Schuler, et. al., 1981). One requirement for such measurements is, of course, 
that the fl<-line can be easily distinguished from random variance in the 
background level. This requirement is frequently not well satisfied and 
7- therefore severely limits applicability as well as the accuracy of the technique 
(Schuler, 1978). 

' The ratio of intensities of Ak-line and background spectrun, x« has been 

shown to be related to the surface wave slope spectral density S(K w ) by the 
equation 

2ir 2 |m| 2 S(K )coth 2 (K d) 

' 1 LJ U7 


where d is water depth, A is illuninated area, and m is the modulation 
transfer function relating Bragg wave modulation to the slope v f the long 
surface wave. (Alpers and Hasselmann, 1978; Plant and Schuler, 1980). The 
inverse dependence of x on A has been observed experimentally (Plant and 
Schuler, 1980; Johnson et al., 1982) and the equation has been shown to 
account well for observed values of x (Schuler, et al., 1982; Johnson and 
Weissman, 1983). Alpers and Hasselman (1978) pointed out that this signal 
to background ratio can be increased by low pass filtering but their best 
estimate of such a ratio was only 1.8 for measurements from satellite 
altitudes. Thus, once again, the background clutter spectrum was the 
predominant factor limiting system performance. 


In this study we have attempted to apply frequency agility techniques 
to dual- frequency scatterometer measurements of the ocean. In this case it 
is necessary to discuss two distinct types of frequency separations. In 
i." addition to the Af separation between transmitted pairs of signals, the 

carrier, or center, frequency of successive pairs of signals will be offset 
from one another by some 6f. The sensor still remains a "dual-frequency" 

-7* system in the sense that the first step in processing is to multiply returns 
corresponding to transmitted frequencies separated by Af. 
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2. EXPERIMENTAL INSTRUMENTATION AND CONDITIONS 


The experimental results presented in this paper were obtained at a 
Naval Research Laboratory field-site located on 36.6 m cliffs overlooking 
the Chesapeake Bay. The dual-channel, coherent radar was built at the Naval 
Research Laboratory and was operated at a frequency of 1.28 GHz. The basic 
operating parameters for this scatterometer system are given in Table 1. 

Both channels of the radar system are fully coherent, i.e., signal phase 
information is preserved. The signal coherency is retained by beating the 
return signal down first to IF, and then to a variable audio offset 
frequency (25 Hz for these experiments) which is set sufficiently high to 
avoid fold-over of any of the spectral components of the backscatter Doppler 
spectrum. Thi3 frequency-agile, dual- transmission channel radar system 
derives its operational flexibility from the nse of an EPROM (Erasable 
Programmable Read-Only Memory) controlled fast-switching (20 us/step) 
frequency synthesizer. The EPROM contains instructions to program the 
synthesizer to create a set of output carrier frequencies [f.] where, i = 1, 
.... N (K N < 98). External counting circuits, controlled by the system 
PRF pufse3, are used to step through the EPROM stored instructions and to 
continually repeat the cycle. The value of N may be manually set to any 
integer value within its range. The individual output frequencies f^ are 
held constant for 100 vs, and then are reprogrammed tc a new value. A 
diagram showing the time dependence of transmitted frequencies corresponding 
to the two received channels is given in Figure 2. These diagrams resemble 
two "staircases" with one offset in time relative to the other by an 
interpulse period. Backscattered returns due to transmissions when the 
channel separation was Af were utilized in the computer processing of the 
data while those due to Af' were suppressed. A synchronization signal was 
derived by using digital counting circuits and a D/A converter. For the 
measurements reported here the difference frequency Af between the two 
L-Band channels (f^, f. + Af) was held constant at a value of 10 MHz. This 
separation allowed us to resonate spatially with water waves of 15 m 
wavelength. This value was chosen because waves of this wavelength often 
are the dominant wave at the Chesapeake Bay site. The values of f^ were 
stepped from 1.230 GHz to 1.328 GHz in 1 MHz steps. The 1 MHz spacing was 
selected a priori using the initial assumption (Pidgeon, 1967) that 
frequency changes equal to or greater than (T; (where T = the transmitted 
pulse width) would be adequate to decorrelate the backscattered returns. The 
work of Voles ( 1966 ) indicates that larger separations may be necessary in 
the presence of long waves but no experimental evidence was to substantiate 
that conclusion. 

The processing of the fully coherent output signal S (t) from the radar 
consists of three steps (1) range gating of the returns, (2) separation of 
the two interleaved frequency channels, and finally (3) processing of the N 
frequency-agile signals which occur during the transmission time of the 
synchronization staircase. 

For the present studies, samples of the backscattered return were taken 
at a two-way delay of 5 us (0.75 km) from the radar site at a grazing angle 
of 3°. The scattering cell dimensions were determined, at this low grazing 
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angle, by the pulse width (lus) and the azimuthal beamwidth of the antenna 
pattern (6° at the two-way 3 dB power point) . These parameters indicate 
that the scattering cell was 150 m radially by 79 m azimuthally. The data 
reported here were taken during an eight hour period on 15 April 1982. 

Wind-generated waves traveled towards shore in water which was 25 m in 
depth; our work indicates a dominant wavelength of 15m. The wind during the 
tests was from the NNE at 15-20 kts as measured by an anemometer at a height 
of 43 m above the sea surface. Tne antennas were pointed into the dominant 
wave direction which was observed to be from the NE. Scattered vhitecaps 
were present. The outputs of the two radar channels as well as the 
synchronization signal were recorded on an analog, four-channel tape 
recorder operating at 30 ips. Subsequent experiments (November, 1983) 
conducted from the CERC Pier Duck, NC have yielded results in agreement with 
the measurements reported on in this paper. 

3. DUAL FREQUENCY DOPPLER SPECTRA USING FREQUENCY -AGILE SIGNALS 

The output data that were recorded consisted of a burst of up to forty- 
nine multiplexed, frequency-agile samples of coherent sea backscatter for 
both radar channels, and a synchronization signal. Our initial processing 
of this multiplexed time signal was both direct and simple. Hie output data 
from the two different frequency channels we-e passed through a laboratory 
multiplier and then through a low-pass filter. The time-constant of the 
filter was set so as to sum the N sequential frequency-agile samples 
together coherently. This smoothed signal was then Fourier transformed by a 
laboratory analyzer (Nicolet Model 444) and a power spectrun was formed. A 
number, generally 32, of these power spectra were then averaged together. 
This mode of processing is fast and effective (i.e., dk-line detectability 
is greatly improved) . Figure 3 illustrates the improvements in the ratio of 
Ak-line to background spectral intensity that can be obtained with twenty 
minutes of sequentially recorded data which had N = 1, 24 and 49 
frequency-agile transmitted signals. It is apparent that the Ak-line 
processed in the conventional way (i.e., with N = 1) is almost unusably 
weak, whereas, the one obtained with N = 49 has been improved by almost 6.5 
dB above the background intensity level. This laboratory or "immediate" 
mode of processing, therefore, produces a significant improvement in signal 
detectability. 

A study was conducted to determine if the dual-frequency resonance line 
detectability improvements derived from the use of frequency-agility would 
be sufficient to obviate the need for additional temporal averaging of 
spectra. If such averaging could be eliminated, then a real-time sensor 
could be developed which would monitor Ak-line strength and line position in 
real-time. Real-time processing would allow ocean surface currents to be 
monitored continuously. The frequency-agile (N = 49) data were played 
through the Nicolet 444 analyzer. This analyzer has the capability of 
allowing selection of any of its 400 frequency bins (via placement of a 
cursor) for study. The output of the selected bin was available at an 
output port of the analyzer and was recorded on an y-t recorder. 

Arbitrarily selecting a detection threshold of 6 dB (line strength above the 
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adjacent background) the data tapes were examined to determine what fraction 
of the time this threshold was exceeded. The time histories of the 
magnitude of the peak and the background are shown in Figures 4(a) and (b) . 

Tne resonance peak exceeded this threshold 16.5% of the time when the 
signal was not frequency-agile (N = 1) and the signal varia: oe was 
approximately equal to the mean. When the "immediate" laboi atory processing 
was applied to the same data using N = 49 frequency-agile carrier 
frequencies, the short-term KiG ssc) signal variance was reduced 
significantly and the detection criterion was exceeded during 95% of the 
length of the time record. 

It is useful to represent this laboratory processing mode 
mathematically so that it can be compared later with the processing 
algorithm used on the computer. The signals coming f-om the two system 
output channels (which are proportional to the coherently detected 
buckscatter^ fields) may be represented as E(t, f R ) anc £(t, f + &f) where 
f is the n 1 ' 0 transmitted carrier frequency of the frequency-agile signal 
set. The low-frequency output of the laboratory multiplier may then be 
represented as 

P n (t) = E (t , f n ) *E * ( t , f n + Af) (2) 

where, f = f + 6f, f =1.23 GHz, and 6f is stepped from 0 to N MHz in 1 MHz 
steps , A? = l8 MHz. ° 

The effect of the low-pass filter which follows the multiplier is to 
convolve P n (t) with a filter weighting function W(t), which is the Fourier- 
transform of the pass-band characteristic of the filter, to form a signal V(t) 
given by 


V(t) 


J P n ( t) W(t - t) dt 
— 00 


or, approximately 


V(t) 


T/2 

J P n (T) W(t - T) dT 
T/2 


n=o 


p„<t> 


where, 

T = [2tf f ]“ 1 , AF f = the filter 10dB cut-off width. 


(3) 
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The signal V(t) is then processed in a laboratory Fourier analyzer (Nicolet 
Model 444). The output of this analyzer is a result of Fourier-transforming the 
input signal V(t), squaring it, and ensemble averaging over a nunber of data 
records to form 


l Ju> ! 2 * \k ? " u> 

• IF 

• L 

iWn 

where, v represents a Fourier transformed quantity and u is radian frequency. 

The notations P n (w) means the Fourier transform of the return signal pair whose 
carrier frequency was N MHz above a base frequency of 1230 MHz. 

It should be noted that the smoothing time T of the filter need not be 
exactly equal to the repetition period T R of the frequency-agile signal set. If 
it is made longer then T_ no new terms appear in Eq. (4) until the extension 
exceeds the temporal decorrelation time (at 1.23 GHz) of the sea itself. 

Gq. (4) shows that the simple processing scheme results in an average 
spectrum which may be represented as a sum of auto- and cross-spectra between 
dual-frequency pairs having carrier frequency offsets tf ranging from 0 to N 
MHz. We expect that returns due to transmitted frequencies separated by more 
than an inverse pulse width (in this case by more than 1 MHz) will be highly 
decorrelated. Thus we would expect .lower background signals for cross- than 
auto-spectra in our sun. Since conventional processing consists of sunming only 
auto-spectra, a net ba ckgr ound reduction 
is to be expected for|v(u)l . 

5. PROPERTIES OF CROSS SPECTRA OF FREQUENCY-AGILE RETURNS 

In order to investigate the validity of the results of the previous section 
in detail, we computed cross-spectra of P n and P, and compared them with the 
power spectrun of P^ itself. To accomplish this, the data described earlier 
were processed using a MINC/LSI-11 computer. The procedure consisted of (1) 
multiplying the multiplexed signals together and, (2) demultiplexing the N 
product signals (Figure 5). The result of this operation was a set of 
individual time-series of sanples of the multiplied fields P n (t) for all 
possible values of carrier frequency offset over the range 0 n < n < N. These 
time-series were then (redundantly) re-sampled and stored as disc files by the 
computer. It should be emphasized that these files contained return signals 
P Q (t) which were due to transmissions whose frequency separation was constant 
(TO MHz) but whose carrier frequency offsets varied in 1 MHz steps. Mean square 


L 


( w) P ^( w) 


(4) 


204 



voltage levels of these digitized signals were constant to within 10? up to 35 
MHz separation. Cross- and auto-spectra were then computed using the digitized 
data. 


In most cases, twenty-five seconds of each of two records corresponding to 
different carrier frequencies were transformed using a uniform time window. Hie 
product of the two transforms was computed and stored. A total of N time 
records were stored on disc, where N = 24 for the first data set and N=36 for 
the second. Therefore the nunber of cross-spectra between records whose carrier 
frequencies were separated by 6f which could be averaged together was N-jf if 
Sf is in MHz. This process was repeated four times so that the total nunber of 
cross-spectra we were able to average together was 4(N-6f). For instance, the 
second data set yielded 144 auto-spectra which were than averaged. In one 
instance, discussed below, we transformed the entire 100 second data record and 
summed various combinations of auto- and cross-spectra. In this single case, 
the second data set was augmented to yield a maximun 6 f of 47. 

Figure 6 gives examples of spectra computed using 25 second records. 

Figure 6(a) is the auto-spectrim computed for the first data set with 48 degrees 
of freedom. Figures 6(b) and (c) show the magnitude and phase of the 
cross-spectrun between digital files corresponding to 6f = 1 0 MHz. These were 
computed with 28 degrees of freedom. Apparent in this figure is the decrease in 
background level relative to dk-line intensity which led to the results of 
Section 3 (note the scale change from (a) to (b)). The phase of the 
cross-spectrun is quite erratic due to low nunber of degrees of freedom used in 
the computation. Note, however, that it is near zero in the frequency range of 
the spectral peak where the coherence between the input signals is relatively 
high. 


Figure 7 shows the result of simulating on the computer the frequency-agile 
processing technique described in Section 3* The resolution is higher in this 
figure than the last since the whole 100 second record was transformed. In the 
figure, N represents the total number of spectra averaged. Note that twice 
this number is not the total degrees of freedom since many records were used 
redundantly. Rather the nunber of degrees of freedom is constant at 96. The 
average of auto~spectra shown in Figure 7(a) correspond to the first term on the 
right of Eq. (4), that is to conventional dual- frequency processing. The 
spectrum obtained by averaging all auto- and cross-spectra from the second data 
set is shown in Figure 7(b). This spectrum corresponds to the frequency-agile 
-.echnique output described by Eq. (4). The improvement in signal-to-background 
ratio for this spectrum compared to that of Figure ft a) is 10.2 dB. This is 
larger than that reported in Section 2 becomes only part of the 20 minute record 
used in Section 2 was used here. Finally, an additional improvement in 
signal-to-background ratio of 2.8 dB is obtained by summing only cross-spectra 
as shown in Figure 7(c) and (d). Thus this method of processing results in a 
total improvement over conventional processing of 13 dB, again confirm.' j that 
signal detectability is better for cross-spectra than auto-spectra. 

5. CONCLUSIONS 

We have demonstrated that processing dual- frequency data in a 
frequency-agile mode can provide significant reduction or, in some cases, 
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complete removal of the clutter background limitations on dual-frequency 
operation. The method consists of forming standard, dual-frequency products of 
signal pairs but then performing the final multiplication or spectral density 
calculation using products formed at different carrier frequencies. 

When Doppler spectra for current measurements are the desired product, we 
have demonstrated that a simple frequency-agile processing scheme (the 
"immediate'' laboratory mode) can produce significant improvement in 
signal-to-background ratio of up to 10 dB. 
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(b) Laboratory processing of the frequency-agile, dual- frequency scatteroaeter data. 
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Figure 6 Auto- and cross-spectra of frequency-agile, dual- frequency 
scatterometer output, (a) Auto-spectrum calculated with 48 
degrees of freedom, (b) Magnitude of the cross-spectrun withAf 
s 10 MHz calculated with 28 degrees of freedom, (c) Phase of 
the same cross-spectrun. 
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Computer processing of frequency-agile data, (a) Auto-spectra 
(standard mode), (b) spectrum resulting from sunning auto- and 
cross-spectra for all possible carrier frequencies, (c) spec trun 
resulting from summing cross-spectra for all carrier frequencies 
greater than zero, and (d) spectrun of (c) on a linear scale. 



